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ABSTRACT 

y  A  generalized  theory  of  acoustic  emission  (AE)  is  developed  on  the  basis  of 
the  theory  of  elastodynaaics  and  dislocation  models.  Aconstic  eaission  sources 
are  represented  as  dislocation  sources  and  include  both  discontinuities  of 
displacement  coaponents  and  tractions.  As  AE  waves  are  observed  at  a  stress  free 
surface.  Green's  functions  in  a  half  space  are  obtained.  Fortran  programs  for 
computing  these  functions  for  nonr-Cauchy  solids  are  uaed  to  calculate  AE 
waveforms  from  a  point  crack  and  moving  cracks.  Their  implication  on  current 
attempts  of  determining  source  characteristics  via  deconvolution  is  discussed.  /■ 


INTRODUCTION 

Detailed  analysis  of  acoustic  emission  (AE)  waveforms  has  been  difficult, 
because  of  the  high  frequency  range  of  AE  signals.  Quantitative  evaluation  of  AE 
signals  as  well  as  theoretical  attempts  to  predict  AE  originating  from  sources 
with  prescribed  characteristics  have  been  made.  The  theory  of  AE  still  faces 
difficult  problems,  the  most  serious  of  which  is  the  absence  of  Green's  functions 
for  relevant  geometries. 

The  theory  of  AE  must  be  able  to  specify  the  nature  of  a  source  starting 
from  a  given  displacement  (or  velocity)  history  at  a  defined  point  of 
observation.  Earlier,  theories  of  dislocations  and  elastodynamics  were  applied 
to  simple  analysis  of  AE  generation  /1-3/,  bnt  only  in  an  infinite  medium.  For 
any  AE  analysis,  this  is  unsatisfactory  as  a  stress-free  surface  must  exist  where 
emissions  are  detected.  AE  in  a  half  space  (in  a  semi-infinite  body)  is  a  good 
representation  of  special  experiments  IA,il .  Pekeris  / 6/  obtained  an  analytical 
solution  of  Green's  function  in  a  half  space  for  a  Cauchy  solid  (Poisson's  ratio 
=  0.2S).  Methods  of  generalized  ray,  normal  modes  and  integral  transforms  have 
been  used  to  obtain  a  limited  number  of  solutions  for  a  plate  /7~9/.  In  most  of 
these  calculations,  an  AE  source  was  represented  by  a  force  impulse  or  force 
couple.  While  this  representation  is  appropriate  in  calibration  experiments  that 
utilize  a  force  step,  the  characterization  of  most  AE  sources  requires 
displacement  steps.  In  an  infinite  medium,  the  spatial  derivatives  of  Green's 
functions  are  used  in  conjunction  with  displacement  functions.  However,  the 
presence  of  a  stress-free  surface  makes  this  practice  untenable.  The  spatial 
derivatives  of  Green's  functions  in  a  half  space  or  of  a  plate  cannot  be  given  in 
an  analytic  form  and  require  elaborate  procedures  evf  -  in  numerical  computations. 

In  the  present  paper,  we  summarize  a  generalized  theory  of  AE  /10,11/  for 
the  representation  of  source  characteristics.  It  is  based  on  the  integral 
formulation  of  elastodynamics  and  the  dislocation  theory.  We  have  employed 
Fortran  programs  for  the  calculation  of  Green’s  functions  in  a  half  space  for  a 
surface  pulse  and  for  a  buried  pulse.  These  computer  programs  can  also  compute 
Green's  functions  of  the  second  kind,  which  are  suitable  for  the  applications  to 
AE  waveform  lation.  Several  dynamic  cases  are  considered. 


We  present  a  generalized  theory  of  AE  using  the  integral  representation  of 
a  solution  in  elastodynaaics  / 11,12/.  Let  D  denote  a  doaain  occnpied  by  a  given 
body  in  the  three  dimensional  space  and  S  denote  its  boundary  surface.  With  the 
assuaption  of  linear  isotropic  and  homogeneous  elastic  body,  an  elastodynaaic 
problea  is  to  solve  the  following  equation  in  D; 

•  e 

t)]  =  (X  +  lOu^  ^  +  l»uit  jj  “  P«i  “  0,  (1) 

where  u^x,  t)  is  a  displaoeaent  field  at  position,  z,  and  tine,  t,  and  the  coaaa 
indicates  a  differentiation  (a,  »  du./dx^).  L[  ]  represents  a  differential 

operator  and  is  used  to  siaplity  expressions,  X  and  |i  are  Laae  constants  and  p  is 
the  aass  density.  Since  the  effect  of  a  body  force  usually  is  not  doainant  in 
elastodynaaic  probleas,  we  oait  the  tera  of  a  body  force  in  equation  1. 


The  solution  of  equation  1  is  subject  to  the  initial  conditions  of  a 
quiescent  past  on  D  +  S  and  boundary  conditions  on  S,  Boundary  S  consists  of  Sj 
and  S2,  where  displaoeaent  g.  is  given  on  S..  and  traction  h.  is  given  on  S.,  as 
follow:  1 


ui(x, t)  =  gi(x,t) 


on  S, 


T[ui(x,  t )  1  *  +  ^ui,  jnj  +  °j,inj*  =  0,1  s2. 

Here  n.  is  the  outward  normal  vector  on  S.  Note  that  T[  ]  is  also  used  as  a 
differential  operator  that  describes  a  relationship  between  a  displaceaent  field 
and  a  traction. 


For  two  arbitrary  displaceaent  fields  u^*,  t)  and  vi(x,  t)  in  doaain  D,  the 
reciprocity  theorea  of  elastodynaaics  states 

/(Llu^x,  t)]*vi(x,  t)  -  LCvi(x,  t)]*ui(x,  t)  )dV 


-  /(Ttu^x,  t)]*vi(x,  t)  -  T[vt(x,  t)]*u^(x,  t))dS, 
S 


(2) 


where  *  aeans  a  convolution  integral  with  respect  to  tiae.  Green's  function  G. 
is  then  defined  as  a  solution  of  the  following  equation:  ^ 


(3) 


LtGij(x,t;x',f)]  =  -  6^  6(x-x')6(t-t* ) . 

Setting  u^  ns  «  solution  of  equation  1  and  to  be  G. .  in  equation  2,  ve  obtain 
the  fol loving  integral  form  as  a  solution  of  equation!. 


•^(x.t)  -/  (Gik(r,  t;*',t')*tk(x',t')  -  Tik(x,  t;x*. t'  )*uk(x',  t' )  )dS. 

8 

where  we  define  Tlu^  -  tk  and  T[Gik3  *  Tik.  Tik  indicates  a  traction 

assoeisted  with  a  displacement  field  of  Green's  function  G,.  and  soaetines  is 
called  Green's  function  of  the  second  kind.  It  is  expressed  as  follows: 

Tij  ‘  ^ij, j*k  +  »*Gik,  jaj  +  “Gij.knj 


Next,  consider  s  domain  containing  a  dislocation.  Ve  assume  homogeneous 
boundary  conditions  on  S  (g.  *  o  on  S.  and  hA  *  0  on  S-)  and  consider  another 
boundary  P  (dislocation  surface)  as  shown  in  Pig.  1.  Ve  apply  equation  4  to 
domain  D  surrounding  boundaries  Sj,  s2  and  P,  and  the  following  equation  results: 


^(x.t)  -/  (G.k(x.t;x',t')*fk(x',t')  +  Tik(x.tfx'.t')*[uk(x',t*)])dF. 
P 


(5) 


Equation  5  represents  any  kind  of  dislocation  sources,  and  provides  the 
generalized  representation  of  AE  source  mechanisms.  Generally  speaking, 
discontinuities  of  displacement  and  of  traction  exist  on  a  dislocation  surface. 
However,  it  is  likely  that  one  or  the  other  has  a  dominant  effect  on  AE 
waveforms.  Therefore,  equation  5  can  be  simplified  to  contain  either  the  first 
or  second  term  in  the  integrand,  or 


Vx,t))  *  /  df  /  Gik(x;x»,t-t')fk(x',t)dF, 


(Sa) 


^(x.t)  =  /  dt'  /  Tik(x;x»,t-t')[uk(x',t')]dF. 

—  F 

Denoting  the  elastic  constant  as  Cwr#,  Tik*Uk3  can  be  expressed  as 


(5b) 


T.t*[uJ  =  C  *G,  [u  ]n  . 
lk  1  kJ  pqrs  ip, ql  rJ  s* 


Unfortunately,  G  it  impossible  to  calculate  analytically  unless  the  medium  is 
infinitely  bounded.  In  cases  of  interest  to  AE  analysis,  it  is  therefore 
imperative  to  use  equation  5b  by  directly  calculating  T..  numerically.  When 
surface  area  F  can  be  regarded  as  infinitesimal  compared  with  domain  D,  the 
surface  integral  in  equation  5  is  evaluated  only  at  a  source  point  x'  and  is 
equal  to  AF.  Depending  on  the  types  of  dislocations,  equations  5a  and  5b  with 
the  initial  conditions  of  a  quiescent  past  can  be  simplified  to  the  following 
convolution  integrals: 


*£(*,1) 


ui(x, t) 


'  AF  /  G..(x;x',t-t')f.(x',t')dt', 

0  * 


AF  /  T  (x;x',t-t')[uk(x',t')]dt*. 
0 


□a 


A  discontinuity  of  the  traction  component,  f.,  in  equation  6  is  equivalent 
to  a  point  force.  Consequently,  equation  6  is  employed  to  analyze  wave  motions 
subject  to  an  applied  force,  such  as  the  breakage  of  a  glass  capillary  and  a 
pencil  lead.  This  equation  is  nsefnl  in  the  transducer  calibration,  bnt  not 
appropriate  for  the  analysis  of  AE  waveforms  dne  to  mierofracturing  events 
(although  it  can  be  nsed  with  added  difficulties).  For  the  latter,  it  is 
preferable  to  nse  equation  7.  Here,  a  discontinuity  of  displacement  component  is 
directly  related  to  the  formation  of  a  crack  or  any  dynamic  movement  in  a 
material.  By  using  equation  5b  or  7,  we  can  account  for  a  moving  dislocation 
that  represents  an  incremental  extension  of  a  precrack  using  the  same  method  as 
the  fault  model  in  the  synthesis  of  seismograms  /13/. 

Considering  the  method  of  AE  observation,  we  need  solutions  of  equations  5 
to  7  at  a  stress-free  surface.  In  order  to  analyze  AE  waveforms  by  equations  6 
and  7,  we  find  that  a  solution  of  equation  3  is  a  Green's  function  in  a  half 
space  which  can  be  substituted  into  the  two  equations.  In  real  experiments,  a 
propagating  medium  is  a  finite  body  so  that  corresponding  Green's  functions 
cannot  be  obtained  easily.  However,  wave  motions  in  a  half  space  are  obviously 
observed  before  reflected  waves  arrive  at  the  observation  point.  Thus,  the 
present  method  can  provide  the  initial  parts  of  AE  waveforms  except  in  very  thin 
plates  and  complex  structures. 

OBRRN'S  CTKCTTONS  IN  A  HALF  SPACE 

The  problem  of  determining  the  elastic  disturbances  resulting  from  a  point 
force  in  a  half  space  is  known  as  Lamb's  problems.  Green's  functions  in  a  half 
space  are  only  available  as  numerical  solutions,  and  these  solutions  cannot 
readily  be  applied  to  problems  of  interest  in  AE  studies.  The  programs  for 
computing  Lamb's  solutions  were  given  elsewhere  /ll/.  By  using  these  programs, 
vertical  surface  motions  of  a  stress-free  surface  due  to  a  step  function  force  on 
the  same  surface  or  that  due  to  a  buried  source  have  been  calculated  1 11/. 


AE  WAVEFORMS  POE  TO  CRACKS 


*.«.  .  A  EsisLSSMSk 

Using  the  program  for  a  buried  source  with  revised  external  functions,  we 
computed  a  Green's  function  of  the  second  kind,  T^,  We  simulated  AE  waveforms 
using  this  solution  and  equation  7.  The  dislocation  model  chosen  for  this  study 
is  the  case  of  a  tensile  crack  parallel  to  the  surface,  or  a  Node  I  crack.  The 
unit  normal  n^  of  the  dislocation  surface  F  is  identical  to  the  x^-axis,  and  a 
displacement  of  the  dislocation  has  only  a  [u^]  component.  From  equation  7,  the 
resulting  displacement  is  expressed,  as  follows: 


•3(x,t)  -  T33<x.t;x't)*(u3(x'.t')J  =  <W33i,i+XG32,2+U  +  2“)G33 , 3} *In3] 


(8) 


The  spatial  derivatives  of  Green's  functions,  such  as  G,j  j  ,  Gj,  *nd  G33  3 
are  computed  separately,  again  revising  external  function^.  In  older  to 
investigate  the  applicability  of  this  method,  we  computed  the  epicenter  response 
of  T.3>  which  is  due  to  s  displacement  discontinuity  of  a  step  function.  The 
result  is  shown  in  Fig.  2. 


function  / 14/,  we  simulated  AE  waveforms: 

i  t  i  A 


By  using  the  following  t 

[»(*>]  -  «o«3(f  *  -  f)  ,lB  (f  *  "  ?)  (° 

Results  are  ahown  in  Fig.  3.  The  rise  time  x  was  assumed  to  be  750  ns.  An 
epicenter  response  (x3  «  2.4  cm)  and  a  response  at  x.  -  6  mm  and  x3  *  2.4  om  are 
ealenlated  to  examine  the  effect  of  a  shift  of  a  source  or  the  location  of  a 


Time  (x  10“®  sec) 


Fig.  4  A  simulated  waveform  due  to  the  moving  or aok  model . giving  [u«] . 

(The  rice  time  l.ps,  a  length  of  the  dieleeation  0.2  mm  and  the 
rupture  velocity. 500  m/s). 

means  of  equation  5a  or  6.  The  crack  moves  0.2  mm  in  the  x^-direction  of  at  a 
uniform  velooity  of  500  m/s.  This  is  represented  by  four  sequential  applications 
of  a  source  displacement  on  the  crack  surface  at  four  points  in  the  x^— direction. 
The  same  source  function  as  in  Fig.  3  was  employed  and  the  rise  time  was  1  |is. 

The  result  is  shown  in  Fig.  4.  In  comparison  to  the  displacement  waveform  due  to 
a  single  crack  given  in  Fig.  3,  we  can  see  that  the  presence  of  a  moving  crack 
broadens  the  peaks  of  the  displacement  response  and  decreases  the  anpl itude  of 
the  S-wave.  The  comparable  experimental  result  of  Wadley  and  Scruhy  /15/  agrees 
with  the  displacement  curve  for  the  moving  crack  quite  well. 


DISCUSSION 

Source  characteristics  of  acoustic  emission  have  been  investigated  as  the 
inverse  problems  using  the  deoouvolution  analyses  / 15/ .  In  these  studies.  Green's 
functions  of  s  point  orsck  were  used  to  deconvolve  AE  waveforms  observed  at 
certain  points  of  observation.  The  effect  of  a  moving  crack  can  be  significant 
on  the  inverse  problem  and  needs  to  be  examined.  As  we  showed  in  the  preceding 
section,  the  moving  crack  broadened  the  AE  waveform.  When  a  faster  crack 
velocity  was  used,  the  AE  waveform  was  closer  to  thst  of  Fig.  3.  Obviously,  one 
must  use  the  dynamic  Green's  function  in  deconvolution  analysis  of  any  observed 
AE  waveforms.  In  the  conventional  deoonvoluton  analysis,  one  is  forced  to  use 
the  Green's  function  of  a  stationary  point  crack  and  it  is  impossible  to  take 
into  account  the  effect  of  a  moving  crack  or  dislocation. 

la  order  to  investigate  the  extent  of  errors  due  to  the  dynamic  nature  of 
cracks  on  the  inverse  problems,  we  performed  deconvolution  analyses  of  simulated 
AE  waveforms.  Simulated  waveforms  were  transformed  into  the  frequency  domain  by 
fast  Fourier  transform  (FFT)  and  were  divided  by  the  corresponding  Fourier 
spectrum  of  Green's  functions  of  the  second  kind  for  s  stationary  point  crack. 
These  deeonvoluted  waveforms  in  the  frequency  domain  were  transformed  into  the 
time  domain  by  inverse  FFT.  Four  examples  are  shown  in  Fig.  5.  These  are 
deconvolved  souree  functions  obtained  from  four  cracks.  One  was  a  stationary 
point  eraek  and  the  other  three  were  fysuie  eracks  with  the  rupture  velocities 
of  2000  m/s,  1000  a's,  and  500  m/s.  Amplitudes  are  normalised  in  those 
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Pig.  s  Source  time  functions  recovered  from  simulated  AE  waveforms  for 

the  moving  cracks.  For  the  cracks  with  the  rupture  velocities  V 
of  2000,  1000,  500  m/s  and  a  stationary  crack. 


The  deconvolved  waveforms  for  the  stationary  crack  and  those  of  two  fast 
rupture  velocities  recovered  the  essential  features  of  the  source  function. 
However,  gradual  decreases  in  displacement  amplitudes  were  observed  beyond  the 
original  rise  time  of  1  ps.  This  decrease  wa  exaggerated  in  the  case  of  the  slow 
crack  with  V  of  500  m/s.  In  this  case,  the  amplitude  became  strongly  negative 
at  3  |is.  This  behavior  apparently  arises  from  the  summations  of  waveforms  from 
different  sources  staggered  in  apace  and  time.  Since  one  has  no  a  priori  * 

knowledge  of  the  source  waveform,  conventional  deconvolution  using  Green's 
function  of  a  stationary  source  may  lead  to  unrealistic  source  functions. 


CONCLUSIONS 

1.  A  generalized  theory  of  AE  is  presented  in  this  paper.  Applying  the 
reciprocity  theorem  of  elsstodynamics  to  a  domain  containing  a  dislocation, 
displacement  fields  due  to  two  components  of  the  AE  source^ function  are  expressed 
by  two  integrals.  One  represents  AE  due  to  an  applied  force  step.  The  other 
represents  AE  due  to  a  discontinuity  of  displacement  components  on  the 
dislocation  surface,  which  corresponds  to  a  crack  or  slip. 

2.  In  order  to  analyze  realistic  conditions  of  AE  detection.  Green's 
functions  of  the  second  kind  in  a  half  space  are  numerically  calculated  by 
Fortran  programs.  Several  representative  cases  are  investigated,  including  a 
stationary  (Mode  I)  crack  and  a  moving  crack. 

3.  Simulated  AE  waveforms  from  dynamio  sources  can  be  easily  obtained 
using  the  present  apprlach.  Displacement  response  from  a  moving  crack  is 
calculated. 

4.  Commonly  used  methods  of  deconvolution  of  sooustic  emission  waveforms 
can  produce  grossly  misleading  conclusions.  While  a  point  crack  response  can  be 
deconvolved  successfully,  the  deconvolution  of  a  moving  crack  responses  may  lead 
to  wrong  source  characteristics.  More  extensive  analysis  of  the  "forward" 
problems  should  be  made  before  attempting  the  "inverse"  problems. 
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